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Abstract. Photons generated by spontaneous parametric down conversion (SPDC) 
are one of the most useful resources in quantum information science. Two of their most 
important characteristics are the purity and the indistinguishability, which determine 
just how useful they are as a resource. We show how these characteristics can both 
be accessed through Hong, Ou and Mandel (HOM) type interferences using a single 
pair source. We also provide simple and intuitive analytical formulas to extract their 
values from the depth of the resulting interference patterns. The validity of these 
expressions is demonstrated by a comparison with experimental results and numerical 
simulations. These results provide an essential tool for both engineering SPDC sources 
and characterizing the quantum states that they emit, which will play an increasingly 
■ important role in developing complex quantum photonic experiments. 

o; 

The ability to engineer and characterize quantum photonic states is of paramount 
importance as the complexity of the experiments in quantum information science 
increases. In particular, it is fundamental to be able to combine multiple independent 
sources and this is usually done via interference, e.g. through a Bell state measurement 
(BSM). At the heart of interference there are two concepts, two characteristics, of the 
photonic states: Purity - how close the real photons are to the theorist's ideal of a single 
mode quantum object, and the Indistinguishability - how alike are the different photons. 
Indistinguishability is a fairly clear concept and its role in interference experiments has 
been widely accepted for some time, in contrast to the purity, which is often ignored. 

The indistinguishability and purity can be increased either by filtering, or by 
engineering the source itself. As we will see, filtering is an option for some applications, 
but it is also a loss factor and potentially a significant one. If one wants efficient, even 
scalable, multi-photon systems, then one needs to go to the source. In this article we will 
focus on spontaneous parametric down conversion (SPDC) as this provides one of the 
most widespread and flexible resources for experiments in quantum information science. 

SPDC sources probabilistically generate two-photon states with diverse features in 
space and frequency and in those degrees of freedom the photons can be anywhere 
from maximally entangled to totally uncorrelated. A common way to characterize 
SPDC generated photons is by using two-photon interferometers. These were first 
introduced by Hong, Ou and Mandel [I] as a tool to measure the relatively narrow 
bandwidth of SPDC photons. The basic idea is to observe the coincidence statistics at 



On the Purity and Indistinguishability of Down-Converted Photons 



2 



the output of a 50/50 beamsplitter as a function of the delay between two photons at 
the beamsplitter input. Due to their bosonic nature the photons may bunch, giving rise 
to the characteristic dip in the coincidences. This interference effect will depend of the 
indistinguishability and purity of the photons and whether they are from the same, or 
independent, sources. Therefore, the question arises as to what information concerning 
the indistinguishability and purity can be extracted from these HOM-like interference 
patterns? 

By focusing on the frequency degree of freedom, in this paper we discuss what 
information can be obtained from the interference patterns of different interferometers. 
We provide simple analytical formulas relating the SPDC parameters with the 
characteristics of the interference patterns. Such formulas can therefore be used to 
engineer photonic states with specific values of indistinguishability and purity. These 
are shown to be in good agreement with both experiment and numerical modeling. 



1. The two-photon state 

In SPDC, pairs of photons are generated by the interaction of a classic pump beam and 
a nonlinear material. Conservation of energy and momentum set the relations between 
frequency and spatial distribution for the pump and the generated photons: signal and 
idler. The correlations between the degrees of freedom can be minimized by a careful 
design of the source [2], or can be erased by projecting the state into a single mode in 
one of the degrees of freedom. In those cases, the two-photon state is a product between 
a spatial and a frequency part. Here, we focus on the frequency part, but an analogous 
analysis is possible for the more general case. 

When considering only the spectral component, the two-photon state generated in 
SPDC is described by the pure state p = |^ r )(^ r |, where: 

|#> = J dnjnM^^(n s )a^n t )\n s ,n t ). (i) 

If we define w n as the frequency of the photons, then Q n = w n — w® gives the deviation 
from the central frequency u>°. All of the information about the frequency distribution of 
the generated pair of photons is described by the normalized mode function $(0,,,^). 
The state of a single photon is calculated by tracing out the other photon from the 
two-photon density matrix p, i.e. p s = Tr [|^ r )(^ r |] i , which can be written as 

Ps = J dM.dtf 9 <i> a (n a ,tt a )\n a )(tt a \, (2) 

where the single-photon mode function is 

Mn s ,n' s ) = J rf^$(fi s ,a t )$*(fi;,a). (3) 

To explicitly calculate the mode function, consider a Gaussian-pulsed laser with 
bandwidth a p , and a crystal with a length L, in that case <&(fi s , fli) is given by 
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$(fi a ,fti) = iVexp 



4^ 



smc 



A k L 



(4) 



x exp [— iAkL] exp 



4(7? 4(7? 



Here AMs a normalization constant, such that |3>(O a , fii)| 2 = 1 and the first exponential 
is due to the spectral distribution of the pump, assuming a Gaussian profile. The 
properties of the material appear through the phase-matching term A& = k p — k s — ki 
and the last exponential represents any possible spectral filtering of the photons, which 
we model as Gaussian functions with bandwidths a s and a^. 

Some approximations are useful when calculating the properties of the mode 
function in (T4J). For instance, by taking the first order terms in the Taylor expansion of 
A*,, it is possible to write the phase-matching term as A^ = (N p — N S )Q S + (N p — N^fli, 
where N p , N s and Aj are the inverse group velocities of the pump, signal and idler. 
Another common approximation is to replace the function sinc(a;) by exp (—ax 2 ) where, 
if a = 0.193, as both functions have the same FWHM. Under these approximations the 
mode function becomes, 



$(fi s ,fii) = N exp 



a„ 2 6 n2 c ~ - 



—nj - -nj - -n s Qi - -(mO s + nfi<) 



(5) 



where we define: 



a = a 2 m 2 + l/a p + 1/a 2 ; (6) 
b = a 2 n 2 + l/a 2 p + 1/a 2 ; 
c = a 2 mn + 1/a 2 . 



The remaining factors m = L(N P — N s ) and n = L(N P — N{), are dependent on the 
crystal length and the difference in the inverse group velocities for the different fields. 

Writing the mode function as a product of exponentials as in (jSJ) makes it possible 
to analytically calculate many quantities related to the two-photon state [2]. It also 
provides a useful and intuitive geometrical picture of what is happening in the frequency 
space of the joint function - the first three terms clearly describe an ellipse with the axes 
defined by a and b and the c component defining a rotation about Q s and fV The last 
term describes the relative phase between Q s and fij. 

In the following sections we derive the functions describing the indistinguishability 
and purity of photons for these sources after passing through two-photon, HOM-type, 
interferometers. 
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Figure 1. A typical HOM interferometer, where a pair of photons generated, in 
this case in a Type-II SPDC process, are first spatially separated by a polarizing 
beam splitter (PBS) before being recombined on a beamsplitter (BS) to interfere. 
The coincidence probability R cc presents an interference pattern when measured as a 
function of the delay time At between the photons arriving at the beamsplitter (BS). 



2. On HOM Interference and Indistinguishability 

The indistinguishability / of two photons quantifies the overlap between their quantum 
states, i.e. for states with no overlap, 1 = 0, while for indistinguishable photons 1 = 1. 
The most natural way to test this, is to try to make the photons bunch, which is exactly 
what Hong, Ou and Mandel did in their seminal paper about two-photon interference 
[1] . In a so-called HOM interferometer, as shown in figure [lj the signature interference 
dip can be observed if the detectors cannot resolve the coherence time of the photons 
[3] . If this holds, the probability R cc of having a coincidence between the outputs of the 
beamsplitter (BS) is given by the integral over the detection times (r a , r b ) of the second 
order coherence function G^(r a ,n) = (pE^\T a )Et\r b )E { b +) (r b )E { a +) (r a )), where E a 
and E b are the field operators at the detectors a and b jU E], that is 



Rcc = J dr a dr b G i2 \r a ,T b ). (7) 

By calculating this probability as a function of the temporal delay At between the 
photons, one can show that the coincidences have the characteristic interference pattern 
given by 

R cc (At) = ~ [1 - I(At)} , (8) 

where 

i(At)= J dn s dni^(n s ,ni)^*{ni,n s )exp[-i{n s -n^At]. (9) 

The spectral indistinguishability between the signal and idler corresponds to the value of 
this function for At = 0, / = 1(0) [6]. Intuitively, the fact that the indistinguishability 
depends on the product of $(fi s , fli) and fl s ) shows that 1(0) quantifies the effect 

of exchanging one photon with the other. 

In practice it is more convenient to extract the indistinguishability from the 
visibility of the interference pattern in ([8]). In this case, the visibility is defined as the 
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difference between the maximum and the minimum values of the interference pattern 
normalized to the maximum value, that is 

_ max[i? ce (At)] - mm[R cc (At)] 

max[i? cc (At)] ' ^ ' 

The maximum coincidence probability is obtained for a large At and is equal to |, while 
the minimum occurs at At = 0, and it is given by |(1 — 1(0)). So, for the interference 
pattern described by (JHJ), the visibility is equal to the indistinguishability: V — I. 



Coincidences (10 s) At' 

80 r ' — ! 




-15 10 

Temporal delay (ps) 



Figure 2. The exact and approximated calculations for the HOM dip compared with 
a dip obtained experimentally. 



If the state of the downconverted photons is described using (jSJ), the integral 
in ([9]) has analytic solutions. In term of the parameters for the SPDC process, the 
indistinguishability between the photons can now be simply given by 



Aab - Ac 2 



1/2 



(a + 6) 2 -4c 2 J ' vll) 
such that the photons are indistinguishable when a = b. The interference pattern 
obtained with the interferometer in figure [1] is a dip described by the function 



R cc (At) = - 

where both the width / and the depth a + b — 2c of the dip are functions of the 
pump beam bandwidth, the filters, and the group velocity mismatch between the 
photons. In non-degenerate configurations, where the group velocity mismatch is non- 
zero, At' = L(N S — Ni)/2 is responsible for shifting the minimum of the dip. 



I exp 



t\2 



t (At - At 

' a + b - 2c 



(12) 
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To validate the approximated solution of ( II 2p . we compare it with the interference 
pattern obtained by solving numerically the integrals of (Q (using the exact mode 
function). We also compared this to experimental data obtained from a Type-II 
configuration where we generated orthogonally polarized photons from a 2 cm PPLN 
(periodically poled lithium niobate) [7] crystal with a 780 nm pump laser beam with a 
0.2 nm bandwidth. The photons are frequency degenerate at 1560 nm and are filtered by 
identical 1.4 nm Gaussian filters. The results are illustrated in figure [2j The estimated 
visibility is 85.5 ± 3.2% for the experimental data, 87.6% for the exact calculation 
and 90.0 % when using the approximations. The difference between the experimental 
data and the theoretical calculations can be explained by a small difference between the 
central frequency of the photons, due to fluctuations of the temperature of the crystal 
during the experiment. This difference is also responsible for the oscillations outside of 
the dip [8]. 



Indistinguishability 




0.0 L 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Signal Filter (nm) 

Figure 3. Indistinguishability as a function of the bandwidth of the filter for the 
signal photon, when the idler is filtered with a 1 nm filter. 

From ([6]), one sees that perfect indistinguishability requires 

\ - \ = a 2 L 2 (N p - Nif - a 2 L 2 (N p - iV s ) 2 , (13) 

and, since in Type-II configuration Ns ^ this condition cannot be satisfied using 
identical filters, unless they are significantly narrower than the natural bandwidth of 
the photons |9j. Therefore, even though it is rather counter- intuitive, to optimize the 
visibility of non-degenerate systems, one needs to use different filters for each photon. If 
we take the same experimental parameters as in figure [2j then figure |3] shows how when 
the bandwidth of the idler filter is 1.0 nm (FWHM) the maximum indistinguishability is 
obtained for a 0.6 nm filter for the signal. For identical filters, with 1.0 nm bandwidth, 
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the indistinguishability is only 94.0%. These values will obviously be different for the 
particular experimental scenario under consideration. 

3. On HOM interference and single photon Purity 

While the indistinguishability is a property of one photon with respect to other, the 
purity is a property of a single photon state. The purity P of a state p, is defined as 
the trace of its density matrix squared P = Tr[p 2 ]. Since it is a second order function 
of the density matrix, a measurement of the purity requires at least two copies of the 
state [10]. One of the most common ways to generate those copies is by using two 
identical crystals, c.f. figure@](a), where two photons, 1 and 2, from independent SPDC 
sources are combined on a 50/50 BS. The interference pattern obtained while varying 
the delay between their arrival times, will have unit visibility as a signature of perfect 
indistinguishability and purity [UJ [12j [13] . 




Figure 4. Variations of the HOM schemes to measure the purity of one photon. The 
two copies of the photon, required for purity measurement, can be generated in two 
different but indistinguishable sources as in (a), or by two consecutive pulses of the 
pump laser as in (b). For this case the temporal distance between the pulses should 
be equal to the pulse-to-pulse time for the laser. 

Although numerous experiments have succeeded in interfering photons coming 
from independent sources, to obtain indistinguishable photons from different crystals 
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is experimentally challenging. Variations on the input mode, the composition of the 
material or its temperature, the output couplings, among many other factors, will 
introduce distinguishability between the photons. We propose a way to overcome this 
problem by generating the two photons in the same crystal by using, for instance, 
different pulses or different time bins. Figure H] (b) shows the principle whereby a 
temporal delay At is set such that photons coming from consecutive pulses arrive 
simultaneously at different inputs of a 50/50 BS. 

Let p\ and p 2 refer to photonic states from consecutive pump pulses 1 and 2, each 
of which is given by ([2]). The input state at the BS is given by the product state 
pi <S> P2- The probability of measuring a coincidence as a function of the delay between 
the photons is 

R cc (At) = X - [1 - P(At)} , (14) 

provided there is one and only one photon in each of two consecutive pulses. The 
function P(At) is given by 



P(At)= I dO 1 dfi 2 1 (fi 1 ,n 2 )0*(fi 2 ,n 1 )exp[-i(n 1 -fi 2 )At]. (15) 

To simplify the notation we use the single-photon mode function given by ([3]). The 
visibility of the interference pattern described by (fT4l) is V — P(0). It is possible to 
show that -P(O) = Tr[pip 2 ], which can be also written as [TH [121 f!3] : 



V = Tr M = ^IMI*^ , (16) 

where ||p|| 2 = Tr[p^ p]. Therefore, the visibility V depends on both the 
indistinguishability between the two photons (\\p 1 — p 2 || 2 ) an d the purity of each of 
them (Tr[p1],Tr[pl]). If it is not possible to distinguish the cases where there is one 
photon in each pulse or, two photons in one of the pulses, the maximum visibility 
decreases to V = P(0)/3. A more general study on the effects of multi-pair emission on 
the interference patterns has previously been detailed elsewhere [HI [151 US E]- 

There are two trivial cases in which these quantities are independent. First, when 
both photons come from a maximally entangled pair, such that Tr[p 2 ] = Tr[p 2 ] = 0, 
whereby the visibility V is a measurement of the photon indistinguishability through the 
norm \\pi — p 2 || 2 . Second, when the photons are totally indistinguishable \\pi — p 2 || = 0, 
and then the visibility is a direct measurement of the purity of the single photon state p\. 
In our scheme, c.f. figure H](b), it is reasonable to assume that an experiment is stable 
from pulse to pulse, such that successive photons are indistinguishable. Therefore, in 
this case, the visibility gives a direct estimation of the purity. 

As before, the simplified mode function of ()5|) allows one to calculate analytically the 
interference pattern described by (jM]) . In the case of indistinguishable input photons, 
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RUM) = ~ 



Pexp[ 



At 2 



;i7) 



where the purity of the input state P = -P(O), can be written as a function of the SPDC 
parameters from Eq. ([6]) as 




P = corresponds to the case where the generated photons are maximally mixed, which 
corresponds to c 2 = ab. The generation of pure photons is achieved for ab » c, which 
is possible by using narrow filters, at the cost of reduced photon flux. Alternatively, one 
can set c = 0, i.e. from flBj): 



-1 = -a 2 L 2 (N p - N S )(N P - Ni), (19) 
a p 



which is only achievable for certain materials [T2J dSl E21 ESI ED] ■ 



Coincidence probability 




indistinguishability 

-20 20 

Temporal delay (ps) 



Figure 5. Interference patterns in the coincidence probability for a two-photon state 
and for two independent photons. In the first case the visibility is equal to the 
indistinguishability of the photons and in the second case it is proportional to their 
purity. We consider the same parameters as used for figure [2l for which the purity is 
49.8%. 
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4. Discussion 

One can see that there is a significant difference between the interference pattern 
obtained in a two-photon interferometer compared to the pattern obtained when two 
independent photons interfere, as described by ( Tl~2l) and ( flTl) and illustrated in figure El 
These two equations are defined in terms of the SPDC parameters and hence provide a 
means to engineer the photonic states with the desired indistinguishability and purity. 
In this paper we have also proposed a novel interferometric technique that allows for 
the purity to be measured using only a single SPDC source. This greatly simplifies 
the experimental complexity, and more importantly, does not rely on any assumption 
that photons from independent SPDC sources are indistinguishable. We expect that 
this simple and intuitive approach to characterizing and engineering SPDC photon 
sources will be of increasing importance for experiments combining multiple sources 
and engineering complex photonic networks. 
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